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In this talk we present a mechanism for leptogenesis which is based on gravity waves produced dur- 
ing inflation. We show that when inflation is driven by a pseudo-scalar field the metric perturbations 
generated during inflation can become birefringent, therefore giving a non-vanishing contribution to 
the gravitational triangle anomaly and sourcing lepton anti-lepton asymmetry. As this asymmetry 
is sourced by the fields which are active during inflation, it is not washed out or diluted by inflation. 
The amount of matter asymmetry generated in our model can be of realistic size for the parameters 
within the range of some inflationary scenarios and grand unified theories. This talk is based on [1] 
which has appeared on the arXiv as hep-th/0403069. 
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I. INTRODUCTION 

One of the puzzles of astroparticle physics and cos- 
mology, which has been around for more than half 
a century and since the existence of antimatter and 
anti-particles established in collider experiments, has 
been to explain, if our universe is mostly made out of 
matter, why and how this has happened during the 
course of the evolution of the universe starting from a 
symmetric soup of matter and antimatter soon after 
the Big Bang. 

The fact that in the parts of the universe visible 
to us there is an excess of matter over antimatter 
has now been backed by the recent determinations 
of the cosmological parameters from the cosmic mi- 
crowave background observations and the WMAP ex- 
periment. Quantitatively this asymmetry is usually 
given through the ratio of excess of baryon density to 
the photon density [2] 



— = (6.5 ±0.4) X 10-1° 



(1) 



where ub = nt — and is the number density of 
photons. This ratio has the nice property that is time 
independent, as the evolution of the ni, and Uj with 
the cosmic Hubble expansion are the same. This is a 
small number, but at the same time it is large enough 
to be a puzzle for models of particle physics. 

About forty years ago Sakharov stated the three 
necessary conditions to generate a matter- antimatter 
asymmetry dynamically from a symmetric initial con- 
ditions [3]: 

i ) our particle physics model should have baryon num- 
ber violating vertices. 
a) CP should be violated. 



Hi) CP and baryon number violating interactions 
should be active at a time when the universe is out 
of thermal equilibrium. 

A baryon excess this large cannot be produced in 
the early universe within the Standard Model (SM) of 
particle physics [4] . This is due to the fact that baryon 
number violating interactions in the Standard Model 
are loop suppressed and the only source of CP viola- 
tion in the hadronic sector is in the Dirac phase of the 
CKM mixing matrix, which is not enough for explain- 
ing the baryon asymmetry observed today. Moreover, 
assuming that the scale of inflation is larger than TeV, 
i.e. the SM is at work after inflation is ended, the 
out-of-equilibrium condition can be created at phase 
transitions or through late decay of massive particles. 
The most attractive choice for a phase transition is 
that associated with electroweak symmetry breaking. 
However, that phase transition is probably not suffi- 
ciently strongly first-order. 

Although what is observed is a baryon asymme- 
try, since the 1980's it has been realized that the 
standard weak interactions contain processes, medi- 
ated by sphalerons {SU(2) instantons), which inter- 
convert baryons and leptons and are thermally acti- 
vated at temperatures greater than 1 TeV. Thus, we 
can also create the baryon asymmetry by creating net 
lepton number at high temperature through out-of- 
equilibrium and CP-asymmetric processes [5, 6]. Sce- 
narios of this type are known as leptogenesis. 

To use the possibility of lepton asymmetry and 
sphalerons we, however, need to fulfill Sakharov con- 
ditions for leptons. Within the usual SM this does 
not solve the issue and we are hence forced to asso- 
ciate the observed baryon asymmetry of the universe 
to physics beyond the SM. 
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Most of the allowed parameter space of the mini- 
mal supersymmetric Standard Model has already been 
excluded and large CP violating phases are strongly 
constrained in supersymmetric models [7] though they 
still could appear in the neutrino Yukawa couplings 
that are used in the Fukugita-Yanagida scenario for 
leptogenesis [6]. Models that can explain the baryon 
excess typically involve exotic nonstandard physics, 
CP violating couplings in the Higgs or supersymmetry 
sectors or in the couplings of the heavy neutral lep- 
tons associated with neutrino mass [8]. In any event, 
there is good reason to seek more effective sources of 
CP- violating out-of-equilibrium physics. 

Here we present a new mechanism for the creation 
of the mat tor- antimatter asymmetry, one associated 
with gravitational fluctuations created during cosmo- 
logical inflation [1] . 



II. OUTLINE OF THE MECHANISM 

Let us first spell out how the three Sakharov con- 
ditions are realized in our model of matter-antimatter 
asymmetry, the gravi-leptogenesis. 

A. Lepton number violation 

The lepton number violation in our model comes 
from triangle anomaly. As it is well-known [9], the 
lepton number current, and hence the total fermion 
number current, has a gravitational anomaly in the 
Standard Model. Explicitly, 



N 
16^2 



-RR 



where 



RR — -^^'^^'^^RafSpaRjS'"^ , 



(2) 

(3a) 
(3b) 



and N = Nl — N^, which is three in the Standard 
Model. The anomaly is a consequence of an imbalance 
between left- and right-handed leptons. In general 
when heavy right-handed neutrinos are also added to 
the Standard Model, as is done in the seesaw mecha- 
nism for explaining the smallness of the neutrino mass, 
(2) will be correct in an effective theory valid below 
a scale /x, of order of the right-handed neutrino mass. 
More concretely N can in general be a function of 
energy. At low energies, below the right-handed neu- 
trino mass scale N = 3. At higher energies, N could 
be anywhere between zero to three, depending on the 
details of the particle physics model invoked. In the 
usual seesaw scenarios with three right handed neu- 
trinos, II can be as large as 10^'*, for energies below /x 



= 3, and above that N = [11, 12]. In any case, 
here we do not restrict ourselves to a specific parti- 
cle physics model, and keep an open mind on larger 
values of /U. 



B. CP violation 

The need for CP violation manifests itself in our 

model through the fact that a non-zero lepton num- 
ber generation can be achieved when (RR) is non- 
vanishing. As we will show explicitly, RR receives 
a contribution with a definite sign from gravitational 
fiuctuations produced during infiation, which is driven 
by a pseudo-scalar field. In other words, CP-violation 
in our model arises from the infiaton field with a CP- 
odd component. Note that during inflation the ex- 
pectation value of the inflaton is non- vanishing and is 
rolling in time. 

Such models of inflation can be naturally achieved 
if the infiaton is a complex modulus field such as one 
finds in supergravity or superstring models. In or- 
der to use these models, however, we need to make 
sure that they have fiat enough potentials required for 
(slow-roll) infiation. The simplest model of this kind 
is when we have a single field inflation and a pseudo- 
scalar (f> as the infiaton, known as natural inflation, 
however it can be incorporated to have multiple ax- 
ions such as in N- inflation models [16]. In fact, these 
models of inflation flt quite nicely into extensions of 
the standard model and in string inspired inflation. 

The imaginary part (f> of this field (which we hence- 
forth call an 'axion') can couple to gravity through an 
interaction 

AL = F{(}))RR , (4) 
where F is odd in (p. Under P and CP 

Terms like (4) would generically appear once we inte- 
grate out heavy fermions axially coupled to </> or as a 
result of the Green-Schwarz mechanism [13]. In the 
Appendix we will explicitly show how a linear F of 
the form 



m = 



(167r2Mpi) 



(5) 



with N depending on the details of string compacti- 
fication, arises from heterotic string theory compacti- 
fied to four dimensions [14]. 



C. Out-of-equilibrium 

Out of equilibrium in our model is achieved noting 
that we apply the interaction in (4) to the dynamics of 
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metric fluctuations during inflation where due to (ex- 
ponential) growth of the background space-time, the 
lepton number production is naturally out of equilib- 
rium. 

We are now ready to compute the amount of the 
lepton number generated during inflation. In our anal- 
ysis, wc assume that wc have a given siicccssful model 
of inflation which involves an axion fleld and do not 
elaborate on the details of the inflationary model. The 
axion driven models of inflation, the natural inflation, 
are preferred from a particle physics perspective. Nat- 
ural inflation gives a small self coupling for the inflaton 
fleld without flne tuning [10]. Throughout this note 
we use notations and conventions of [15], in particular 
Mpi = 2.44 X 10^^ GeV is the reduced Planck mass. 



III. GRAVITY WAVE EVOLUTION 

To begin, we must compute the production of grav- 
itational waves during inflation under the influence 
of the coupling (4). The action which describes the 
gravity waves is hence 



a = -M^^v'-det g R + F{4>)RR. 



(6) 



In general metric perturbations about an FRW uni- 
verse can be parameterized as 



ds'^ = -{1 + 2tp)dt^ + uj.dtdx^ 



(7) 



where ■)/', Wi and hij respectively parameterize the 
scalar, vector, and tensor fluctuations of the metric. It 
is straightforward to show that the scalar and vector 
perturbations do not contribute to and so we ig- 
nore these fluctuations in the following discussion. We 
can also flx a gauge so that the tensor fluctuations are 
parameterized by the two physical transverse trace- 
less elements of hij. For such physical gravity waves 
which are moving in the z direction, the metric takes 
the form 

ds'^ = -df +a^{t)[{l- h+)dx'^ 

+ {1 + h+)dy'^ + 2hxdxdy + dz"^] (8) 

where a{t) = e^* during inflation and h+, h X are 
functions of t, z. 

To see the CP violation more explicitly, it is conve- 
nient to use a helicity basis 



1 

V2' 



Hl = —7={h+ 



ihx) , Hr = -^{h+ + ihx) ■ (9) 



Here /il and hfi are complex conjugate scalar fields. 
To be very explicit, the negative frequency part of /iz, 
is the conjugate of the positive frequency part of hn, 
and both are built from wavefunctions for left-handed 
gravitons. 



A. The equations of motion 

Plugging (8) into the action (6), up to second order 
in Hl and Hr, we obtain 



: -{hLDhn 
16iF{4>) 



^hRDhL) 

— h h 
dz"^ ^ dtdz ^ 



—h -^h 

dz^ ^ dtdz 



dt^ "^dtdz 



^^'^ ^dt^^dtd-z^'^-dt^'^dtd-z 



hR 



(10) 



where 



dt"^ dt a? dz'^ 

As it is explicitly seen from (10), if ft-L and Hr 
have the same dispersion relation, RR vanishes. Thus, 
nonzero RR requires "cosmological birefringence" 
during inflation. Such an effect is induced by the 
addition of (4) to the gravitational equations. Lue, 
Wang, and Kamionkowski (LWK) [18] and Alexander 
and Martin [19] have studied the effects of such an in- 
teraction in generating observable parity-violation in 
the cosmic microwave background. In a future work, 
this coupling will be related to an observable to de- 
tect birefringence in binary systems for the LISA and 
Advanced LIGO gravitational wave detectors [17]. 

It is straightforward to obtain equations of motion 
for h^ and hR-. 

UhL = -'^i^f^L , ^hR = +2i^h'ji , (11) 



where 



~ 4(F"^2 ^ 2F'H4))/Mli 



(12) 



dots denote time derivatives, and primes denote differ- 
entiation of F with respect to (/>. To obtain the above 
equations we have used the fact that the inflaton fleld 
is only a function of time t. In the second line of the 
expression for 0, assuming the slow-roll inflation, we 
have dropped the terms proportional to cj) (explicitly 
we have dropped iF'cJy/Mpi). 

Note that with a constant 0, (4) is the Gauss- 
Bonnet term and being a total divergence (in four 
dimensions), cannot affect the equations of motion; 
thus, all terms in Q involve derivatives of (p. These 
equations should be compared to those for evolution 
in flat space given by LWK [18]. The new term pro- 
portional to Hcj) leads to a substantial enhancement 
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in the size of G. With this simphfication, and the 
approximate form (5), 



27r2 \Mpi 



(13) 



where e = ^{^Y / {H M-pxf' is the slow-roll parameter 
of inflation [15]. 



B. Gravitational birefringence 

To see gravitational birefringence we need to solve 
the equations of motion explicitly. Let us focus on the 

evolution oi Hl and, more specifically, on its positive 
frequency component. It is convenient to introduce 
conformal time 



1 



1 



-m 



Ha H 



(Note that conformal time rj runs in the opposite di- 
rection from t.) The evolution equation for Hl then 
becomes 

—.hL-2-—hL--r^hL = -2iQ——hL. (15) 

If we ignore for the moment and let h^ ^ e**^^, 
this becomes the equation of a spherical Bessel func- 
tion: 



rf2 Id 
—hL-2- — . 
ar}^ T] ar] 

for which the positive frequency solution is 
/i+(fc,ry) = e+*'=(''+^)(l-ifcr?) . 



(16) 



(17) 



We now look for solutions to (15) with Hl e'*^^. 
To do this, let 



hL = e''' ■ (-ifcr?)e'=®''ff(r?) 



(18) 



where g{r]) is a Coulomb wave function, we then have 
d^ 



d^n'''' 



9 = 0. (19) 



This is the equation of a Schrodinger particle with 
^ = 1 in a weak Coulomb potential. When 6 = 0, 
the Coulomb term vanishes and we find the spheri- 
cal Bessel function (17). For h^, the Coulomb term 
is repulsive; for Hr, with the opposite sign of the Q 
term, the Coulomb potential is attractive. This leads 
to attenuation oi Hl and amplification of ha in the 
early universe. This is the anticipated cosmological 
birefringence which was also discussed by LWK [18]. 
(For a more detailed treatment see [19].) 

As we will see generation of the matter asymme- 
try is dominated by modes at short distances (sub- 
horizon modes) and at early times. This corresponds 



to the limit kr] ^ 1. In this region, we can ignore 
the potential terms in (19) and take the solution to 
be approximately a plane wave. More explicitly, 

5(77) = exp[zA:(l - 6^)^/2^(1 + a(r?))] , (20) 

where 0(77) ~ log 77/77. 

IV. THE GREEN'S FUNCTION 

We would now like to use (18) to compute the ex- 
pectation value of RR in an inflationary space-time. 
It turns out that it will be dominated by the sub- 
horizon, quantum part of the gravity-wave evolution. 
Hence, to compute (RR) we only need the two point 
(Green's) function {hLhn): 



(14) G{x,t;x',t') = {hL{x,t)hR{x' ,t')) 



I 



'^'^^g./c.(.-.')G,(,7,77') . ^^^^ 



For k parallel to z, the Fourier component Gk satisfies 
(15) with a delta- function source 



r j2 



^-2(^; + ^0)i; + fc^ G.(77 , 77' 



dr] 



-5(77-77') 



For 6 = 0, the solution of this equation is 

m+{k,v)h^{-k,v') V<ri' 

miik,r])h+{-k,7]') rj'<r] 
where 

H V 1 



(22) 



M 



PI 



2fc3 



(23) 
(24) 



and is the complex conjugate of (17), and h~^, 
are the corresponding solutions of the Hr equation. 
For 6 = 0, these solutions are the same as for h^, but 
the structure of (23) will be preserved when we go to 
the case 6^0. The leading effect of 6 is to introduce 
the exponential dependence from (18), 



Gk 



-k&ri +k&ri' 



GkO 



(25) 



for both rj > rj' and r] < rj' . The prefactor is modified 
in order 6^, and the wavefunctions acquire additional 
corrections that are subleading for krj ^ 1. Neither 
of these effects will be important for our result. 

The Green's function (25) can now be used to con- 
tract hL and Hr to evaluate the quantum expectation 
value of RR. The result is 



{RR) 



(26) 
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where we have picked up only the leading behavior for 
kr] > 1. 

We would like to emphasize that our expression for 
{RR) is nonzero because of the effect of inflation in 
producing a CP asymmetry out of equilibrium. The 
original quantum state for the inflaton might have had 
nonzero amplitude for a range of values of and might 
even have been CP-invariant. However, inflation col- 
lapses the wavefunction onto a particular value of (j) 
that is caught up in the local expansion of the uni- 
verse. This value gives us a classical background that 
is CP-asymmetric. 

The above result and computations seems to be cru- 
cially dependent on the form of the Green's function 
or the vaciiiuii state wc; have used. To resolve the 
possible ambiguity in this regards, one may perform 
the above computation using a different method, the 
fermion level crossing, e.g. following [20]. This com- 
putation confirms the above results [21]. 



V. COMPUTING THE RATIO OF LEPTON 
TO PHOTON DENSITIES 

To complete our leptogenesis model, there are two 
remaining steps. First, we need to compute the net 
lepton anti-lepton asymmetry generated. This, in our 
model, should be computed during inflation. We then 
need to have a reheating model through which we can 
compute the temperature and hence the entropy of 
the Universe after the inflation ended. 



A. Lepton number density 

Wc are now ready to evaluate the lepton density 
that arises through the gravitational anomaly (2). In- 
serting (26) into (2) and integrating over the time pe- 
riod of inflation, we obtain 



n ■ 



Jo ''J(^ 



d^k 1 SH^k^rfQ 



)3 167r2 



Ml, 



(27) 



where n is the lepton number density. The integral 
over k runs over all of momentum space, up to the 
scale /X at which our effective Lagrangian description 
breaks down. The dominant contribution comes not 
from the modes which have left the horizon by the 
end of inflation (the super-horizon modes), k/H < 
1, but rather from very short distances compared to 
these scales. The integral over rj is dominated at large 
values of r;, toward the end of inflation (c/. (14)). The 
integral represents a compromise between two effects 
of inflation, first, to blow up distances and thus carry 
us to larger physical momenta and, second, to dilute 
the generated lepton number through expansion. It is 
now clear that the dominant contribution to the right- 
hand side comes from krj ^ 1, as we had anticipated. 



We can now perform the integrals to find the lepton 
number density produced by the end of inflation 



1 



n = 



727r4 \MpiJ 



(28) 



Let us now analyze each factor in n: 

• The factor (H/Mpi)^ is the magnitude of the grav- 
ity wave power spectrum. Wc should stress that the 
usual gravity wave power spectrum is coming from the 
super-horizon modes, while the main contribution to 
n has come from the sub-horizon modes. 

• The factor is a measure of the effective CP viola- 
tion caused by birefringent gravity waves. 

• The factor is the inverse horizon size at inflation; 
this gives the density n appropriate units. 

• Finally, the factor (/i/i?)^ gives the enhancement 
over one's flrst guess due to our use of strongly 
quantum, short distance fluctuations to generate RR, 
rather than the super-horizon modes which effectively 
behave classically. 



B. Photon number density 

To understand the significance of the lepton num- 
ber density (28), we should compare it to the entropy 
density of the Universe just after reheating, or to the 
photon number density. (Recall that almost all of the 
entropy of the Universe is generated during the re- 
heating time and is carried in the massless degrees 
of freedom, i.e. photons.) In order this we need a 
reheating model. 

To illustrate the physical relevance of the result of 
our model, we take the simplest (and at the same 
time very naive) reheating model, the instant reheat- 
ing model. That is, we assume that all of the energy of 
the inflationary phase has been converted to the heat 
of a gas of massless particles instantaneously. The 
reheating converts the energy density of the inflaton 
fleld 6 



to radiation with 



(29) 



(30) 



and to the entropy density s, s = 271^ g^,T^/ 45, 
where T is the reheating temperature, is the ef- 
fective number of massless degrees of freedom and 
s = 1.85* Uj [15]. (Here we have assumed that the 
evolution of the Universe after the reheating era has 
been adiabatic). This gives 

= 1.2Sg-^/\HMpif/\ 

Recalling that the ratio of the present baryon number 

to the lepton number originally generated in leptogen- 
esis is approximately ns/riL = 4/11 [5], in our model 
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we obtain 



— = 4.05 X 10" 



5 3/4 

y* 



Mpi) 



e 



(- 



(31) 



If we are less naive, we might follow the dilution of n 

and p with the expansion of the universe to the end of 
reheating. The final result is the same (sec, however, 
[23] for a comment on this point). With the adiabatic 
expansion assumption, (31) can be compared directly 
to the present value of ns/n^ given in (1). 



C. Numerical results 

We should now estimate the numerical value of 
UB/n^ in our model to see if we have a viable lepto- 
gcnesis model. In our final result (31) we have five di- 
mensionless parameters, 5*, H/Mpi, iJ,/Mpi and the 
slow-roll parameter e and K (the latter two appear 
through the CP violation parameter G). Inserting the 
expression for Q we have 



riB 



~ 2.9 X 10 



— y 



-1/2 



M 



PI 



(32) 

Within the usual supersymmetric particle physics 
models 51* ~ 100 is a reasonable choice. The WMAP 

data, through the density perturbation ratio 5p/ p (for 
a single field inflation) leads to an upper bound on 
H/Mpi ratio as [22] 



H 



< 10" 



e < 0.01, 



(33) 



ox H< 10^^ GeV. 

The factor Jvf in (5) is inferred from the string the- 
ory compactification and is proportional to the square 
of the four dimensional Mpi to ten dimensional (fun- 
damental) Planck mass [14] (see also the Appendix). 
Within string theory. 



10 



K~ 10^ - 10 



could be a reasonable range. Therefore, assuming that 
the scale of inflation H saturates its current bound 
H/Mpi ~ 10-^, 



e ~ 10"^ - 10°. 



(34) 



The physically viable range for the parameter fi de- 
pends on the dc;tails of the underlying particle physics 
model on which our gravi-leptogenesis is based. For 
example within the Standard Model + three heavy 
right handed neutrinos (and the seesaw mechanism) p 
could be of order of the right handed neutrino mass, 
in which case p, < lO^^GeF. If we do not restrict 
ourselves to the seesaw mechanism, p can be larger. 
The upper bound on p is then only coming from the 
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range of vaUdity of our effective field theory analy- 
sis. For example, we can accommodate /z's as large as 
the string scale where quantum gravity effects become 
important. Therefore, depending on the details of the 
model 

IQ^^GeV <p< 10^^ - lO^GeF 



or 



10-6 < 



M, 



< 10-2 - 10- 



pi 



(35) 



is a reasonable range for p. 

In sum, assuming that scale of inflation H saturates 
its current bound H ~ lO^^Ge^, we have 



TlB 



10+^ e 



Mpi 



(36) 



Noting (34) and (35), it is readily seen that there 
is big parameter space for obtaining the observed 
baryon asymmetry of the Universe within our gravi- 
leptogenesis model. 



VI. SUMMARY AND OUTLOOK 

Here we have constructed a leptogenesis model in 
which the lepton asymmetry is a result of gravita- 
tional chiral anomaly. Indeed our model is a "mod- 
ule" , which uses minimal ingredients and could be fit- 
ted into a successful inflation model in which the in- 
flaton field(s) has a pseudo-scalar component. Prom 
the particle physics side, it is again like a module and 
many different particle physics models could be in- 
voked. Here we did not discuss the details of neither 
the inflationary nor the underlying particle physics 
models. It is of course very important to explicitly flt 
our module in viable physical models. 

As our construction is a leptogenesis model, we 
need a mechanism to convert the lepton asymmetry to 
baryon asymmetry. With the Standard Model this is 
usually performed via the thermally activated electro- 
weak sphalerons, which in our model should be active 
after inflation. This happens for reheat temperature 
Tr > 1 TeV. Within our naive "instant reheating" 
model that translates to H > 10~^ eV, which is not a 
serious constraint on the inflationary models used. 

According to recent WMAP observations, the scalar 
metric perturbations generated during inflation have a 
size that gives density fluctuations with 5p/p ^ IQ-^. 
On the other hand, noting (1), one has 



Ub 



(37) 



The above could be a simple numerical accident or 
there might be some underlying deeper physics. In 
any case, our model which invokes gravity waves as 
the source for the baryon asymmetry and operates 
during inflation, could be an interesting framework to 
uncover the possible deeper physics behind (37). 
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Appendix A: Realization of the model within 
string theory [14] 

To see how this term appears in string theory, con- 
sider e.g. the Heterotic SUGRA action: 

S =Mfo j di°a;v/d^t^(3l + '^{d^f 

where 

{^^{A) and ^^{co) are the gauge and gravitational 
Chern-Simons three-forms, respectively. Explicitly, 

2 

i^3{oj) = Tr{oj Adoj + -u Aoj Aoj) 



In particular note that 

*{dn3{uj))=RR. 

Upon compactification to four dimensions, the 
term leads to 

S ~ M^o j d^y j d'^x e-'^{dB + a'n3{Lo)f 
= Mfo— / d^x[{dBf + 2a'{*dB) A ^3(0;) + • • • ] 

= ^ / d'^x{{d4>f - 2a'4>RR) 
9s J 

where the pseudoscalar </> (our inflaton field) is dual 
to the two form B2 and the 4d Planck length, Mpi is 
defined as 



M2 = MfoVe. 



And hence, finally: 



9s \MioJ 
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